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We use first-principles density-functional total energy and 
polarization calculations to calculate the piezoelectric tensor 
at zero temperature for both cubic and simple tetragonal or- 
dered supercells of Pb3GeTe4. The largest piezoelectric coef- 
ficient for the tetragonal configuration is enhanced by a factor 
of about three with respect to that of the cubic configuration. 
This can be attributed to both the larger strain-induced mo- 
tion of cations relative to anions and higher Born effective 
charges in the tetragonal case. A normal mode decomposition 
shows that both cation ordering and local relaxation weaken 
the ferroelectric instability, enhancing piezoelectricity. 



Piezoelectric materials are of great technological im- 
portance. The highest known piezoelectric coefficients 
are found in mixed compounds, such as Pb(Zri_3;Ti3;)03 
near its morphotrophic phase boundary at x « 0.45. This 
association has been reinforced by the recent discovery of 
giant piezoelectricity in single crystals of the relaxor fer- 
roelectric systems Pb(Ai/3Nb2/3)03-PbTi03, (A = Zn, 
Mg). Q Understanding of the physical origin of enhance- 
ment of piezoelectricity in mixed ferroelectrics would be 
of great theoretical interest and could also point to new 
ways to tune this property. 

Ab initio calculations have proved successful in relating 
properties of ferroelectrics to phenomena on the atomic 
level. Structural parameters Mjm, dielectric constants, ef- 
fective charges and phonon dispersion relations Mpl , and 
structural phase transition temperatures M have been 
calculated from first principles. Such methods can nat- 
urally be extended to piezoelectricity. Piezoelectric co- 
efficients have been calculated from first principles for 
various materials |3"|9[|, where good agreement with ex- 
perimental results is obtained. For mixed systems, ab 
initio calculations can be performed on individual or- 
dered supercells and the results analyzed to understand 
the dependence on configuration and the resulting en- 
semble average. In this Letter, we consider the mixed 
ferroelectric Pbo.75Geo.25Te. Calculation of the piezo- 
electric response shows significant configuration depen- 
dence, with enhancement of piezoelectricity in the cell 
with lower symmetry. We interpret and explain the re- 
sults in terms of chemical bonding and long range dipolar 
interactions. 

Pbi_2,Gea;Te undergoes a transition from a paraelec- 
tric cubic rocksalt phase at high temperature to a ferro- 



electric rhombohedral phase at low temperature, dom- 
inated by Ge off-centering fin], for all compositions 
X > 0.005. |ll|. In order to study the dilute Ge limit 
while maintaining small supercells, we choose x = 0.25. 
The configurations which we investigate are shown in 
Figure ^. The structures have identical compositions 
Pb3GeTe4, but one is simple cubic (cP8) and one is sim- 
ple tetragonal (tP8). 
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FIG. 1. Pb3GeTe4 cP8 and tP8 configurations. 

First-principles density-functional calculations were 
performed using the program CASTEP 2.1 [O, with 
a plane wave basis set in the local density approxima- 
tion. Bachelet-Hamann-Schliiter pseudopotentials |1^ ] 
were used for all atoms. For the cP8 configuration, we 
used a 4 X 4 X 4 Monkhorst-Pack set of k points; for the 
tPS configuration, we used the same k point grid, folded 
into the tetragonal Brillouin zone. In our study of the 
structural phase transition in the Pb3GeTe4-cP8 config- 
uration |14| , we give further details on the ab initio calcu- 
lations. Fixing the atoms in rocksalt positions and mini- 
mizing total energies with respect to strain, we obtained 
a = 6.275 A for cP8 and a = 4.438 A; c = 12.548 A for 
tP8. These are taken as the high symmetry reference 
structures for the two configurations. Next, we deter- 
mined the ground state structures. For a series of fixed 
global strains, as many as 6 independent internal coordi- 
nates in the cP8 configuration and 16 internal coordinates 
in the tP8 configuration were relaxed via conjugate gra- 
dients minimization. In each relaxation, we broke the ini- 
tial symmetry by displacing each Ge atom by 0.3 A in the 
(111) direction. A quadratic fit of total energies as a func- 
tion of strain was used to find the equilibrium strain. The 
cP8 ground state obtained is rhombohedral, space group 



R3m, a = 6.315 A, a = 89.47°, and the tP8 ground state 
is monoclinic, space group Pm, a = 4.497 A, 6 = 4.409 A, 
c = 12.575 A, /3 = 89.31°. The internal coordinates are 
given in Table |, with the center of mass fixed in each 
case to that of a strained rocksalt structure with Ge on 
the origin. 

The method of King-Smith and Vanderbilt (KSV) |l|] 
was then used to calculate polarizations. This method 
calculates the electronic contribution to the polarization, 
to which the ionic contribution is added to obtain the to- 
tal polarization P. The component of polarization along 
lattice vector a.^ is expressed in terms of a set of wave- 
function calculations along straight-line paths in the Bril- 
louin zone, where reciprocal lattice vector be is the path 
vector, J the number of k points in the path, and kj^ the 
component of k orthogonal to b^ along the path. 

In Voigt notation, the components of the (third rank) 
piezoelectric strain tensor are given by Cap = dPa/dep, 
where Cq. is a component of the (second rank) strain ten- 
sor. We used finite differences to determine the piezo- 
electric tensor components: Ca/j ~ APa/Ae/j. For both 
the cPS and tP8 configurations, we calculated P for the 
ground states and for a set of fully relaxed strained cells 
with one strain component e^ increased by 0.005. Strains 
were measured with respect to the axes shown in Fig- 
ure ^ Strain magnitudes were measured with respect 
to the cP8 reference structure (a = 6.275 A) in both 
cases. The KSV calculations were performed using the 
self-consistent potentials obtained from total energy cal- 
culations using 64 k points in the Brillouin Zone. We 
used J = 32 for cP8, J = 16 for tP8 along bs, and 
J = 48 for tP8 along bi and b2. In each case, we ob- 
tained the cross section average by averaging the values 
for kj_ on a 4 X 4 Monkhorst-Pack grid. These param- 
eters were sufficient to converge the piezoelectric tensor 
components to within 0.05 C /m? . 

The resultant piezoelectric tensor components and 
their decomposition into clamped ion and internal strain 
contributions ^ arc shown in Table ||. Generally, the 
components in the tP8 configuration are larger in mag- 
nitude than in the cP8 one. The clamped ion terms are 
remarkably similar; the differences are almost entirely 
due to internal strain. We will focus on 633 because it 
is the largest piezoelectric component in each case and 
shows the greatest change in going from the cP8 to the 
tP8 configuration. 

We first look at the configuration dependence of the 
internal strain term essjint by analyzing the contributions 
of individual ions: 
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where Z* are the Born effective charge tensors and 
dui/de^ is the movement of the relaxed ion position from 
the clamped ion position with strain (in fixed center of 



mass coordinates). For each configuration, we deter- 
mined Z* from the changes in polarization resulting from 
displacing each atom in turn 0.01 A along each axis direc- 
tion, and dui/des by comparing the relaxed coordinates 
of the ions in the two cells differing only in £3. Since 
the displacement derivatives dui/de^ are dominated by 
motion in the z direction, wc make the following approx- 
imation: 



633 lint « 633 lint = y 2_^ ^i,zz-Q^ 



unit cell 

E 42 lint (2) 



The calculated values for du^z/de^ and Z*^^ are given 
for each ion in Table |lll[ 

For both the cP8 and tP8 configurations, the values of 
633 lint obtained from Eq. || agree with e33|int to better 
than 6%. The values of Table III can thus be used to 



associate the enhancement of piezoelectricity in tP8 with 
the response of the individual ions. It can be seen that 
the increased ionic motion under strain is the dominant 
factor, with additional enhancement due to the larger 
average magnitude of the Born effective charges. Both 
the magnitude of the piezoelectric response and the con- 
figuration dependence is greatest in the chains along z 
containing the Ge ions. 

The theory of resonant p-bonding ||l^,|l^ in IV- VI com- 
pounds gives insight into this trend. In the symmetric 
rocksalt structure, the covalent bonds are composed pri- 
marily of p orbitals, with one electron per bond. In each 
anion-cation- anion segment, there is thus the tendency 
for the cation to prefer an off-center position, forming 
a short cation-anion bond on one side into which the 
electron from the long bond is partially transferred. The 
off-centering is therefore associated with a large Z*. This 
tendency is counterbalanced by short range repulsion, the 
strength of which is roughly determined by the ionic ra- 
dius of the cation. In pure PbTe, short-range repulsion 
dominates and the low temperature structure is rocksalt 
||l7| , while in GcTc, the off-centering instability domi- 
nates and the low temperature structure is rhombohc- 
dral. In Pbo. 935 Geo. 015, the low temperature structure is 
rhombohedral at low pressure and cubic at high pressure 
||l8| , showing that the relative importance of the short- 
range repulsion increases with pressure. 

Equivalently, a decrease in the Te-cation-Te distance 
suppresses off-centering. Figure shows the z-component 
of the polarization of pure GeTe and PbTe as a function 
of lattice parameter, fit to the results of KSV calcula- 
tions in which internal relaxation is allowed but the lat- 
tice is held cubic. A ferroelectric instability is predicted 
for PbTe at sufficient negative pressure. In the mixed 
system, we define local lattice parameters along z to be 
equal to the Te-cation-Te distance, and qualitatively as- 
sociate the local piezoelectric response with the slope of 
the Pz vs. a curve for the corresponding pure cation- Te 
system. In the cP8 configuration, there is a single local 



lattice parameter, the value of which favors off-centering 
and piezoelectric response for the Te-Ge-Te chain but not 
for the Te-Pb-Te chains. In the tP8 configuration, there 
is a different type of chain, with alternating Ge and Pb. 
In this chain, relaxation of the Te ions along z leads to 
two effective local lattice constants: a short Te-Ge-Te and 
a large Te-Pb-Te distance marked in Figure 0. The lat- 
ter Te-Pb-Te distance is considerably larger than in pure 
PbTe, favoring off-centering and piezoelectric response in 
the Te-Pb-Te segment as well as the Te-Ge-Te segment, 



as seen in Table III 
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FIG. 2. Polarization vs. lattice parameter in PbTe and 
GeTe and first-principles LDA "local lattice parameters" 
along z for PbTe, GeTe and mixed Pb3GeTe4 configurations 

Another informative way to analyze the piezoelectric 
response is to decompose the ionic displacements into the 
contributions of the zone-center normal modes of the sys- 
tem. For the decomposition of 633 |mt, we need only the 
polar modes of the high symmetry reference structures 
involving ionic motion along z. For the tP8 systems, two 
reference structures were considered. One, tP8„„, has 
all ions sitting on ideal rocksalt coordinates, the other, 
tPSrei, is derived from tP8„„ via a symmetry preserv- 
ing relaxation. In this relaxation, the Te atoms in the 
z — ±0.25 plane move =f0.143z A and the Pb atoms in 
the z — ±0.25 plane move =p0.029i A. The z polar modes 
were determined from first principles frozen phonon cal- 
culations, as in Ref. |14|, with the results given in Ta- 
ble [V. For each reference structure, one mode has imag- 
inary frequency, indicating instability to a symmetry- 
breaking distortion. 

The approximate piezoelectric response e^^\int for a 
given ground state can be decomposed into normal mode 
contributions: 



e-SMint 



V 2^ ^^' 563 ' 



where mode effective charges are given hy Z = 
^unit CO ^^^^^ {duiz/dS,^) and d^^/de^ is the change in 
the amplitude of mode /i with strain. All modes are 
normalized so that the sum of the squared atomic dis- 
placements is 1 A^ per 8-atom cell. 

The results are shown in Table IV. In each case, the 
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largest mode contribution is from the unstable mode. 
In fact, the unstable tP8r-ei mode accounts for 97% of 
e^^lmt for tP8. We find that the relative Pb(3) motion 
in the unstable mode goes from 0.16 for tP8„„ to 0.49 
for tPSrei, in agreement with the local viewpoint that a 
sufficient increase in the Te-Pb-Te distance will favor Pb 
off-centering and piezoelectric response. 

We investigated the unstable mode contribution in 
greater depth by constructing a simple "double well" 
parametrization of the energy as a function of the am- 
plitude of unstable mode distortion and strain 63. When 
all other parameters are equal, decreasing the strength 
of the quadratic instability increases the piezoelectric re- 
sponse associated with the unstable mode. This trend is 
clear in Table IV. The weakening of the tPStm instability 
relative to the cP8 instability is the result of the change 
in cationic configuration and is due largely to less fa- 
vorable dipole-dipole interactions. Symmetry-preserving 
relaxation further weakens the instability and leads to an 
even greater unstable mode response. 

As v^ -^ 0^ for an unstable mode, the associated 
piezoelectric response will diverge. In stoichiometric 
ferroelectrics, this is essentially what happens at the 
paraelectric-ferroelectric transition, with a correspond- 
ing divergence in the piezoelectric coefficients on the 
low-symmetry side. In substitutionally disordered com- 
pounds containing one stable and one unstable endmem- 
ber, such as Pbi_a;Gea;Te, there will be a variety of local 
environments and therefore a range in the degree of local 
(in)stabilities. We expect in such systems that enhanced 
piezoelectricity can be observed for a range of tempera- 
tures, pressures and applied fields where the density of 
marginally unstable optical modes is significant. Even 
at fixed composition, the type of ordering present has a 
dramatic effect on the local stability and could be used 
to tune the piezoelectricity, {e.g. through control of pro- 
cessing or in epitaxial thin films). Quantitative mod- 
eling of the temperature-dependent response associated 
with the unstable modes of particular configuration, us- 
ing the effective Hamiltonian approach |19[| , is currently 
in progress. 

In conclusion, we have compared the piezoelectricity in 
two Pb3GeTe4 configurations and found strong enhance- 
ment of 633 in the tetragonal configuration. This effect is 
intrinsic; it will occur even in a single grain, single phase 
sample. The interpretation of the calculated results in 
terms of local polar instabilities and the long-range in- 
teractions between them suggests further opportunities 
for identification of solid solutions with large piezoelec- 
tric response. 
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TABLE I. Ab initio ground state structures of Pb3GeTe4 
cP8 and tP8 configurations. Lattice constants given in text. 
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Atom 
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Pb(l) 


0.0014 


0.5037 


0.5037 


Pb(l) 


0.5120 


0.5 


0.2485 


Pb(2) 


0.5037 


0.0014 


0.5037 


Pb(2) 


0.4961 


0.5 


0.7516 


Pb(3) 


0.5037 


0.5037 


0.0014 


Pb(3) 


0.0054 





0.5065 


Ge 


0.0308 


0.0308 


0.0308 


Ge 


0.0576 





0.0062 


Te(l) 


0.9985 


0.9985 


0.4760 


Te(l) 


0.9985 





0.2349 


Te(2) 


0.9985 


0.4760 


0.9985 


Te(2) 


0.9849 





0.7563 


Te(3) 


0.4760 


0.9985 


0.9985 


Te(3) 


0.4661 


0.5 


0.9980 


Te(4) 


0.4955 


0.4955 


0.4955 


Te(4) 


0.4958 


0.5 


0.4966 



TABLE II. Comparison of piezoelectric tensors for the 
ground states of two Pb3GeTe4 configurations (in C/m ). 
Components in parentheses are equal to other components 
via symmetry. 



Component 


Total 


Clamped ion 


Internal strain || 




cP8 


tP8 


cP8 


tP8 


cP8 


tP8 


en 


1.8 


2.5 


-0.5 


-0.5 


2.3 


3.0 


ei2 


-0.6 


-1.2 


0.0 


0.0 


-0.6 


-1.2 


ei3 


(-0.6) 


-0.9 


(0.0) 


0.0 


(-0.6) 


-0.9 


ei4 


-0.2 


-0.6 


0.0 


0.0 


-0.2 


-0.6 


ei5 


1.0 


1.2 


0.6 


0.5 


0.4 


0.7 


ei6 


(1.0) 


1.1 


(0.6) 


0.5 


(0.4) 


0.6 


631 


(-0.6) 


-0.3 


(0.0) 


0.0 


(-0.6) 


-0.3 


£33 


(1.8) 


5.1 


(-0.5) 


-0.5 


2.3 


5.6 


634 


(1.0) 


2.2 


(0.6) 


0.7 


(0.4) 


1.5 


636 


(-0.2) 


-0.5 


(0.0) 


0.0 


(-0.2) 


-0.5 



TABLE III. Contributions of individual ions to esalint in 
Pb3GeTe4 cP8 and tP8 configurations. Atoms located in 



chains along z containing Ge are italicized. 

C/m^; ^ in A. 
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cP8 






tP8 




Atom 


zu 




633 int 


Atom 


Zi^zz 




633 jint 


Pb(l) 
Pb(2) 
Pb(3) 

Ge 
Te(l) 
Te(2) 
Te(3) 
Te(4) 


5.4 
5.4 
5.9 
3.1 
-4.0 
-5.0 
-5.0 
-5.8 


0.72 
0.72 
0.91 
0.79 
-3.08 
-0.37 
-0.37 
-0.46 


0.25 
0.25 
0.34 
0.16 
0.78 
0.12 
0.12 
0.17 


Pb(l) 
Pb(2) 
Pb(3) 

Ge 
Te(l) 
Te(2) 
Te(3) 
Te(4) 


5.9 
5.8 
5.6 
5.7 
-5.9 
-5.7 
-5.5 
-5.9 


0.32 
0.91 
3.08 
2.60 
-3.11 
-2.80 
-0.83 
-1.74 


0.12 
0.34 
1.11 
0.95 
1.18 
1.03 
0.29 
0.66 



TABLE IV. Polar mode contribution to eaajint in 
Pb3GeTe4-cP8 and tP8. cP8 and tP8 modes have symme- 



try labels Pfs and Pi, respectively. Z^ are in eA; C33I 
C/m? and v^ in cm^^. 
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633 1!::! 


I^M 
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dijdes, 


- |(M) 


!^M 


zl 


a^M/aea 


e^^ 


v^, 


X 


diJde-s 


103 i 


5.9 


2.93 


1.09 


64 i 


11.5 


6.28 


4.64 


40 i 


13.6 


6.32 


5.53 


31 


1.6 


1.09 


0.11 


42 


2.9 


-0.14 


-0.03 


42 


5.2 


0.64 


0.22 


46 


4.2 


2.44 


0.65 


85 


6.6 


3.30 


1.41 


89 


3.1 


-0.62 


-0.12 


82 


10.6 


0.52 


0.34 


102 


3.6 


-1.08 


-0.25 


95 


4.2 


0.46 


0.12 


105 


1.0 


-0.23 


-0.01 


117 


1.3 


-1.20 


-0.10 


122 


4.0 


-0.28 


-0.07 



